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I. INTRODUCTION 

The physics of electronic localization is one of the most interesting issues in the context of disordered quantum 
systems [1]. The localization/delocalization transition has been intensively studied by following the field theoretical 
approach of Wegner [2], and the more formal developments by McKane and Stone [3]. Universal aspects of this 
transition and related topics like transport phenomena in disordered media can also be studied in the framework of 
random matrix theory (for a review see, for instance, ref. [4]). 

In an integer quantum Hall system, electron localization plays a fundamental role [5] , namely, the effect of disorder 
is that of localizing the states of energy E away from the Landau energy En = (n + l/2)?icjc- Near this value, the 
localization length £,{E) diverges like £,{E) (x {E — En)~'^ . Although this behavior is consistent with experiments [6] 
and with several numerical [7] and renormalization group calculations [8], to the best of our knowledge, there is no 
simple analytical method to evaluate this scaling [9] nor it is clear why, in the thermodynamic limit, just one state 
for each Landau band remains extended. 

The electronic properties are encoded in the system's correlators; for instance, the density of states (localized or 
not) can be obtained from the retarded Green's function, 

p(£:) = i lim Im{G{x,x;E-ir])), (1.1) 

where (...) denotes the ensemble average over the random variable. On the other hand, the analytical structure of 
the density of localized states pi{E) is more complex since it involves the mixed product of retarded and advanced 
Green's functions, 

pi{E) = - [ dx \im v{G{x,y;E ~ irj) Gix,y;E -\- irj)). (1.2) 

In fact, any quantity encoding some information about localization or transport properties is related with this mixing. 
For this reason, we concentrate our efforts in studying the representation of the ensemble average for these important 
quantities. 

There are essentially two methods to calculate averages: the replica [10] and the supersymmetric method [11]. 
The first one is based on the introduction of a number of copies to represent each type of Green's function. In 
order to compute the average density of states one set of n field copies is needed. In this case, a generalized n 
component Landau-Ginsburg model is obtained [12], where the limit n — *■ must be implemented. For the density 
of localized states two sets of copies are used, containing n+ and n_ field copies, respectively. The disorder average, 
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after integrating out all massive modes, leads to a generalized nonlinear sigma model with symmetry group 0(n+, n_) 
[2, 3], describing the critical properties near the mobility edge. 

In the supersymmetric method, the Green's function denominator is exponentiated by means of Grassman variables, 
thus avoiding dangerous continuations in the number of copies. In order to compute the density of states a superfield 
is needed, this also applies to other quantities that can be computed from correlators just involving a retarded 
prescription. This method has been used by Brezin, Gross and Itzykson [15] to study the density of states in Hall 
systems. They showed that the reason behind the exact expression previously calculated by Wegner [2] is an exact 
boson-fermion symmetry in the Green's function representation, when the projection to the first Landau level is 
considered. However, the problem is not so simple when one tries to extend these ideas to understand localization. 
As occurs in the replica method, a doubling is used to represent the density of localized states, that is, a superfield is 
introduced for each prescription in eq. (1.2). As a consequence, the boson-fermion symmetry is lost, thus precluding 
the possibility of obtaining exact results (see the discussion in ref. [15]). 

Therefore, it is interesting to look for alternative representations for mixed correlators, thus opening the possibility 
to advance the understanding of disordered systems. 

The aim of this paper is to propose a new representation for mixed correlators based on the introduction of essentially 
just one superfield. Our method is general and can be applied to any noninteracting disordered system. It is based 
on the introduction of a single operator, quadratic in the system's Hamiltonian, that permits to encode the nontrivial 
information about retarded-advanced mixing by using a single prescription. Then, the basis of the method is in fact 
quite general and could be useful outside the supersymmetric framework. 

Our method will be tested in the simpler context of random matrix theory where the supersymmetric method has 
been introduced by the pioneering work of Efetov [13]. We will also show how the boson-fermion symmetry can be 
restored when considering mixed correlators in Hall systems. 

In section II, wc outline the general single superfield method for mixed correlators. In sections III and IV we 
compute the corresponding averages for the Gaussian Unitary Ensemble, while in section V we study level spacing 
correlators and compare with well known universal asymptotic behavior. In section VI we briefly discuss quantum 
Hall (QH) systems. Finally, section VII is devoted to present our conclusions and discuss possible perspectives. 



II. SINGLE SUPERFIELD/SUPERVECTOR REPRESENTATION 

We start by considering the operator, 

0{£)=H-£, (2.1) 
where H is the system's Hamiltonian. We will use the notation, 

Oi+ = 0{Ei + irj) and O2- = 0(^2 - «??), (2.2) 

where rj > and Ei, E2 are two real energies. Now, we will define the operator Om, which is quadratic in H and 
mixes the retarded and advanced prescriptions, 

Om = O1+ O2- = O2- O1+ (2.3) 
= {H-Ef-{r/2 + ir]f, (2.4) 

where E = {Ei + E2)/2 and r = Ei — E2- As the factors in Om commute, it is easy to verify, 

Om [Or+ - 02_i] = O2- - O1+ = (r + 2i7^), (2.5) 
therefore, the inverse of the operator Om is, 

(for rj ^ 0, these inverses are well defined). In other words, the three types of associated Green's functions G{x, y) = 
{x\0~^\y) are related by, 

Gm{x, y; E, r) = ^^-^ [G{x, y; E^ - ir^) - G{x, y; E2 + iv)]- (2.7) 

Our strategy is based on first using a single superfield to compute four-point correlators associated with the 0,„ 
operator. Then, from Wick's theorem and eq. (2.7), we can relate this correlator to the sum of three terms 
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involving only retarded, only advanced or mixed products of two-point Green's functions, respectively. The first 
two terms are in general simpler to study and well known in many different disordered systems. Then, we see that 
the four-point correlators for Om operators encode important information about mixed products. Moreover, as the 
retarded-retarded and advanced-advanced terms also admit a single superfield representation, the ensemble average 
for the mixed prescriptions could be performed on independent single superfield representations. 

In particular, in the random matrix case, this procedure will be tested by calculating finite N correlation functions 
and comparing with well known large N universal results for the Gaussian Unitary Ensemble (GUE), given by the 
probability distribution, 

P{H)dH = J\fl[dHul[dH,j dH,j exp{~-jN tr{H^)), (2.8) 

i=l i<j 

where H is an. N x N hermitian matrix. 
Taking traces in eq. (2.6) we get, 

tr{0-') = [tr(Or+) - tr{0^-)]/{r + 2iry). (2.9) 

Now we can multiply eq. (2.6) twice, take the trace, multiply eq. (2.9) twice and then sum together the two expressions 
to obtain, 

tri0^l)tr{0^1) + tr{0-,l0^1) = 

= l/2[tr{0^l)tr{0^l) + tr{0^l0^l) + tr{0^})tr{0^}) + tr{0^10^1)] 

-(r + 2z,7)V2[tr(0„i)tr(0„i)+ ir(O-iO-i)]. (2.10) 

When rj goes to zero, the first member in cq. (2.10) is a mixed retarded and advanced four-point correlator at 
energies Ei and £2, respectively. 

On the other hand, in the second member, we have the sum of three four-point correlators for the operators O1+, 
O2- and Om, respectively. For each correlator we can introduce a single supervector representation of the type, 

tr{0-^)tr{0-^) + tr{0-^0-^) = 

(A/2)2(-l)^(^+i)/27r-^ / d(pdipd^pdi^{ip-lpfeM-S), (2.11) 



5=^(^0<p + V0V), (2.12) 

where Lp (tp) is a bosonic (fermionic) A'^-component complex vector. The parameter A must be chosen for the integrals 
to be well defined. 

Then, the idea is using three independent single supervector representations of the type given in eq. (2.11) to 
compute the average for the mixed correlator. This is in contrast with the usual supervector doubling to represent a 
typical term in the first member of eq. (2.10), 

(K+\ij02-\ki = (l/2)^7r~^^ j d(fi+ dAf^ dtp+ dtp^ dip- dxp_ dtp- dtp- (p+iTp_^_jip-k'ip-i 

X exp ^-^{lp_^_Ol+(p+ +'^fJ_^_Ol+'^lJ+) + ^{lp_02-(p- +tp_02--ip-) 

(2.13) 



III. AVERAGES FOR THE FOUR-POINT CORRELATORS OF Om OPERATORS 

When averaging over the imitary ensemble in the second member of eq. (2.10), the more involved part of the 
calculation comes from the Om correlator, as it is associated with an exponent which depends quadratically on H. 
Then, in this section we will concentrate on this part, providing a detailed derivation of the correlator, 

{tr{0-')tr{0-')+tr{0-'0-')) = 

(A/2)2(-l)^(^+i)/27r-^ J d^ dip dtp # {if ■ lpf{eM-Sm)), (3.1) 



5m = ^^OmV'+^VO^V- (3.2) 

Using eq. (2.4), and considering the matrix K = (p ^ip — tp (E) tp, we can also write, 

(exp(-S'„)) = ATexp Q[(r/2 + irjf - E"^] trK^ J dH exp(- trHSH + XEtrHK), 
where S = \^NI + ^K. Completing squares, we have, 

{eM-Sm))=Me^^[^[{Tl2 + inf-E'']trK+'^E''trKS-^K^ (2/7iV)^ J(5), (3.3) 
where we have defined, 

I{S) = I dH exp(- trHSH) , S = I + -\-K. (3.4) 

In order to evaluate I{S), we will consider the general case where S = I + g{K), substituting g{x) = at the end 
of the calculation. 

As the measure dH is invariant under unitary transformations H ^ U HU^ , U G U (N) , we have, I{USU^) = I{S); 
and taking the matrix U that diagonalizes S: USU^ = D, S = WDU we can write, 

I{S) = I{D) = J dHexp ^-Y,diHl - ^(d, + dj)HijHi}j 

= 7r'^l[di'/^l[{di + dj)-^=n'^2^exp j -1/2 ^In^ + d,) j , 

where di = 1 + gi are the diagonal elements of D, written in terms of the eigenvalues gi of the matrix g{K), so we 
obtain, 

I{S) = (V2)*2f exp (|-l/2i:C (^) j > C(^) =ln(l+a=). 
Now, if we expand ({x) = J2n=i o-nX^, the exponent in I{S) can be written in terms of U{N) invariants. 



(3.5) 



'■,3 " 



= r M^'M."-'). (3.6) 

In the appendix, we show that the trace of a given function f{K), with K ~ Lp ®Tp — jp ® tl), is given by, 

trf{K) = Nfo + /(a) - /(/3) + "^'^"^ ~ {^^^ M, (3.7) 

a — p 

a = (fi -Tp , (3 = ip ■ ip , n = ip- ij). 
Using this formula in eq. (3.6), we can resume the series in the index n, 

Ec(^) = 

ag{a)) + 2NC {g{a)/2) + CigiP)) - 27VC(g(/3)/2) - 2C(.9(a)/2 + g{P)/2) + 
+ ^ i ^'^^'^^"^^ + 2iVC(ff(a)/2) - 2C(ff(a)/2 + g{P)/2)] + 

■ ^'^ [agm-2NC{g{f3)/2)-2C{g{a)/2 + g{l3)/2)], (3.8) 



a-P df3 
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where we have used g(0) = and C(0) = 0. 

Now, from eqs. (3.5), (3.8), and also considering the formula (3.7) in the first factor of eq. (3.3), we finally obtain. 



(exp(-5rr»)) = Gi{a,/3)+G2{a,p)ijJl, 



(3.9) 



where, 



'yN + Xa 77V + A/3 



{2 + X0/JN) 



N 



(2 + A(a + 0)hN) 



and 



(2 + XahN)^ (1 + Aa/77V)V2(i + A/3/77V)V2 ' 
(X'^E'^ (2 + A(a + /3)/77V) 



(3.10) 



A2 



\ 2'yN (1 + Aa/77V)2(1 + Xp/-fN)^ ^ -f^N{2 + Xa/jN){2 + Xp/-fN) 
X^{/3-a) 



+ 



273jv3(i + Xa/jN){l + Xp/jN){2 + X{a + P)/jN) 



(3.11) 



Note that this expression is correctly normalized as for ip = 0, ip = wc have, Sm = 0,a = 0,/3 = and /U = 0. 
Then, in this case, we see from eq. (3.10) that the first and second members in eq. (3.9) are equal to one. 
Now, we can substitute eq. (3.9) in the correlator (3.1) to obtain, 

{tr{0-')tr{0-')+tr{0-'0-')) = 

(A/2)2(-l)^(^+i)/27r-^y" d^dipd'4jcbpa^[Gi{a,/3)+G2ia,P)fi-p]. 

(3.12) 

We will still transform this expression to a function of a and f3 only. In order to do so, we recall that the fermionic 
integral picks-up the term of the integrand having the form Z)^j^^i'!/'2V'2 • • • V'jvV'JV- Now, in l3 = ijj ■ ip, the ■^j's and 
V'j's come in pairs V'iV'ii then, the crossed terms {i ^ j) in 

Hji = if ■ ip Tp ■ '4) = ^ ^iTpitpiiJi + ipiTpjiJilpj, 
i ijij 

cannot contribute to the fermionic integral, that is, 

j dil)dil)G2{a,p)nJi j dipdtljG2{a,P)Yipi(Pitjjiil)i 

= jjYl'^^'^i j ##G2(a,/3)t/' • = jj^-'f j dipdil)G2{a,f3)'4} ■ i), 



(3.13) 



where we have also used the symmetry among the ■^j's. Putting this information together, we arrive at, 

(ir(o-i)MO™^) + MO™0-i)) = 

(A/2)2(-l)^(^+i)/27r-^y d(^c^^##a2[Gi(a,/3) - ^G2(a,/3)]. 



(3.14) 



Now, the integrals over ip, can be easily expressed in terms of a single variable integral over a. Using the solid angle 
in M dimensions, £Im = 27r^/^/r(M/2), and the fermionic integral. 



J di;dijF{P) = (-l)J^(^+i 



)/2 



dl3 



N ' 



0=0 



we finally obtain. 



{tr{0-')triO-') + tr{0-'0-')) 



(A/2) 



(TV. 



)i r 

l)!io 



da 



1 d 



3JV 



aj3 



Gi{a,l3)-^G2{a,l3) 



(3.15) 



/3=0 
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IV. ANALYTICAL DETERMINATION FOR THE CORRELATORS OF Om OPERATORS 



Here, we will study the analytical properties of the four-point Om correlators. Recalling that is a small parameter 

1). Wc will 
i. Precisely, 



(wc will take the limit ry 0), wc can drop if' in the exponent ■|(r/2 + irifa, in cqs. (3.10) and (3.11). Wc will 



suppose r = Ey — Ei> Q\ then, for the integral representation (3.15) be well defined, we must take A 
this choice produces an associated factor in the bosonic integral. 



expiwa exp(— ?7ra/2) 



r 

w = — , 



(4.1) 



which leads to a convergent integral. 

We underline here that, as usual, the factor TV in the exponent of the probability measure in eq. (2.8) assures that 
the energy band attain a finite width when large A'^ matrices are considered. Therefore, the natural scale for the level 
spacing is 1/A''. For this reason, we will define the variables. 



s 



(4.2) 



Using eqs. (3.10), (3.11) and (3.15), working around E = {Ei + i?2)/2 = 0, and considering the change of variables 
a ct/ijJ, (3 —ifi/uj, the small rj {r} « r) contribution of the last term in eq. (2.10) to the mixed retarded and 
advanced correlator results. 



C{X) = -Y{triO-')triO-') + tr{0-'0-')) 



Jo 



da a^~^^ exp(ia) 



1 d 



N 



N\dl3 



N 



exp(— /3) X 



/3=0 



(# + 6)^-1 (a + hf 1 (# + br {a Jr)ib " «) 



(X+a)Af+i (afo)V2 4x(^+a) 



AT 



(a6)3/2 



(4.3) 



where a = ^ 



ia, b ■ 



^ + /3. 



JV 



We see that this expression is a sum of terms containing products of the form /!^"^(5^"\ with n = 0, 1 and fj,, v 
semi-integers, where. 



■ 7Z 



da ■ 



a 



N+l( X 



(2^+ia)^+i-" 



and the Q's are polynomials. 



1 d 



N 



N\d(3- 



N 



13=0 



5)^-^+"6'*exp(-/3). 



(4.4) 



(4.5) 



Now, the integrand in cq. (4.4) presents singularities on the upper complex a-planc, on the imaginary axis. There 
is an isolated singularity when a = 2ix/N, and there is a continuum of singularities for a = iy, y > xl^ ■ As no 
singularities are found on the first quadrant, and the a integral on the arc at infinity goes to zero, we can rotate 
the integration path until the positive imaginary axis is approached anti-clockwisc. Recalling that, from cq. (4.1), lo 
contains a small positive imaginary part that regularizes the integral, according to eq. (4.2) we have also to consider, 
X ^ X + Then, we can use the equivalent representation. 



ii"\x) = i(-if^' 



y^'^'^i^ - y + ie)" e:icp{-y) 
(2^-y + ie)^+i-" 



We can split this integral into two parts. 



/(") = J(")+M"), 



(4.6) 



(4.7) 



where the integral in the J's runs from to x/-^i while the integral in the K's runs from x/-^ to oo. A straightforward 
calculation leads to, 



4"HX)=7(-1)'^+'(X/A^) 



n+u+l f 

Jo 



dv 



?;^+i(l- w)"exp(-^v) 



(2-^) 



N+l-n 
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On the other hand, we can write, 

Note that the factor e^^'^ is pure imaginary, as the ly's arc semi-integer. Then, using {2x/N — y + ie)~^ = P{2x/N — 
y)~^ — i-K5{2x/N — y), we see that the real part of the K^s can be easily computed, as this is associated with the 5 
term. 



1 rl^~^ 

SRiirW(^) = e^'^("-"+5)77^- 



y^+'(y-^)^exp(-y), (4.8) 



N 



{N-n)\ dy^-"^ 
and changing variables, y = 2x/N + j3, 

^Kl^\x) = e-(^-"+^) Pi")(x), (4.9) 

where we have introduced the functions, 

2^ _L M^ + l^f 



1 jAf- 

J^i"Hx)=77rexp(-2x/A^). 



{N-n)\ 

and b has been defined after eq. (4.3). 



(^+ 6)^+^6- exp(-/3), (4.10) 

/3=0 -'^ 



V. CORRELATIONS FOR LEVEL SPACING 

In this section, in order to test our single superfield method, we will first obtain an approximate relationship between 
the Om-correlators and the DOS-DOS correlators, which becomes exact in the N ^ oo limit. Then, we will compare 
the finite N O^-correlators, as N becomes large, with the well known universal results for DOS-DOS correlations. 

The density of states operator, at energy Ei , is given by, 

i 

(cf. eq. (2.2)). Then, the DOS-DOS correlator can be written as, 

p{E,,E2) = {p{Ei)p{E2)) - {p{Ei)){p{E2)) = 

^-L_^(^tri0^l)tr{0^^)) - {tri0^l)){tr{0^1))) 

~2^^ ((*KOr+) tr{0^l)) - {trO^DitrO^D) • (5.2) 

As we are interested in analyzing the new single superfield representation for mixed correlators, we will simplify the 
standard part of the calculation as much as possible. In particular, in the N ^ co limit, the following factorization 
of retarded-retarded correlators is verified (see for example refs. [4], [17]), 

{triO^DtriO^D) = {trO^DitrO^l). (5.3) 

In this manner, up to irrelevant terms that scale to zero for large values of A'', we can consider the approximation, 

p{E^,E2)~ 

^-i^JR (-rV2( tr(O-i) tr(O-i) + tr(O-iO-i))) + (5.4) 



R_ 

N2 N2 



(5.5) 
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where we have used cqs. (2.10), (5.3), and r ^ 0, — > 0. 
In addition, we have, 

(MOr+Or-^)) = {Y.{E,-E,-ir^)-^) = ^^{tr{0^l)), 

i 

{triO^lO^D) = ±-{tr{0^})), (5.6) 
while, using eq. (2.9) for r ^ 0, the last term in eq. (5.5) can be separated in only retarded and only advanced parts, 

{tr{0-l0^1)} = litriO-D) - litriO^D). (5.7) 

Note that, using (5.2), (2.10), (5.6) and (5.7), the density-density correlator can be generally written in terms of 
four-point Om correlators, up to only retarded and only advanced terms which can be studied by following standard 
procedures. 

In particular, in the Random Matrix case, we quote the well known large N results, 

1 tr{0-l) = -1e, + zTTpiE,) , 1 tr{0^1) = -1e^ - iirpiE^), 



7 /4 ^2.4 



Then, using eqs. (5.6) and (5.7), we see that the three terms in the second line of eq. (5.5), when multiplied by the 
factor, are irrelevant. In fact, the second line tends to zero faster than 1/A'', as from eqs. (5.6) and (5.7), 



j^^iHO^lO^D) « -^{E, - E,) = -1 

thus, implying that a, 1/N factor multiplying the sum in the parenthesis of this line would already give a vanishing 
large N result. 

With regard to the first line in eq. (5.5), let us recall that, as explained after eq. (4.1), the appropriate scaling to 
study level spacing is ;• = s/N, with s fixed. This is the natural variable wc have used in the mixed correlator in cq. 
(4.3) (x = 7S^/8), which we are studying at Ei + E2 = 0, that is, Ei = s/{2N), E2 = —s/{2N). As a consequence, 
in the first line of eq. (5.5), we can use, {l/N){tr{OY+)) « i7rp(0) and (l/iV)(ir(0^i)) w -inp{0), that is, 

2^ « ^{liiMO)r + li-iMO)r-iiMm-iMO)) 

« -ipm'. (5.9) 
Then, using this information in eq. (5.4), we find. 



^ p{Er)p{E2) ~ ^ (p(0))2 

1 



■3?(-rV2(tr(0^i)tr(0-i) + tr(O-iO-i))) 



27r2(p(0))2iV2' 

where we have used /5(0) = -^/tt (cf. eq. (5.8)) and eq. (4.3). 

In this manner, we see that the large N four-point correlator of Om operators in eq. (5.10) gives the usual 
{p{Er)p{E2))/{p{Ei)p{E2)) correlator. 
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In order to test our single supcrfield representation for mixed correlators, we can compare the large N behavior for 
the last line in eq. (5.10) and the well known (N oo) universal behavior for the left member in that equation, 

(note that = nr/A, where A = l/{p{0)N) is the mean level spacing). 

As we have seen in section IV, C(x) is a sum of terms containing products of the form lu^'^Q^T^ , with real coefficients. 

The polynomials Q^"'' are real (cf. eq. (4.5)), so that 5ft C(x) depends on the real part of the J's, which have been 

splitted into two terms in eq. (4.7), by using the ji"^ and kI"^ functions. It can be seen that the J's contribution to 
5ft C{x) is suppressed for large N. On the other hand, the real part of the 's has been evaluated in closed form in 
eq. (4.9). Putting all this information together, we arrive at, 

~ 27 I 3/2 ^-1/2 ^1/2 -^-1/2 ^^3/2 4^-^-1/2 '^1/2"'" 

"^4^^1/2 "^-1/2 4jyf^-3/2 "^1/2 4JY^l/2^-3/2 4_^2 ^-3/2 "^-1/2 4_^2 -f^- 1/2 ^^-3/2 



(5.12) 



In figure 1, we display, in the range s G [0, 12], the N oo universal law in eq. (5.11) and our = 10 and = 70 
correlator {2^N^)~^ (C(x)) in eq. (5.12). We see that the correlators approach the correct asymptotic form as N 
increases for all values of the scaled variable s. For small values of s the asymptotic law is well approximated even 
with small values of N. As N increases, larger values of s progressively fit the — > oo behavior. 

To show this effect in more detail we present in figure 2 a set of figures in the range s G [2, 12] where we display 
a detailed view of the universal law and the convergence of our finite N correlator, by considering N = 10, 30, 50 
and 70. We can see that for larger values of N more oscillations are accommodated, clearly showing the onset of 
the asymptotic universal law. Wc; also note in this figure, that in a region close to the origin the convergence is 
very fast. The parameter 7, that characterizes the GUE, has been taken equal to one, however, the same matching 
has been also verified when 7 is changed. In this set of figures, we have also displayed the corresponding finite N 
DOS-DOS correlators obtained in terms of orthogonal polynomials, as given in Mehta's book [18]. Of course, the 
fitting of these polynomials with the universal law is also improved as N increases. Although both curves converge 
to the asymptotic universal law, the finite N results for {2'yN'^)~^ 5i(C(x)) do not coincide with the finite A'' exact 
DOS-DOS correlators. This is expected; as we have seen in section V, both quantities differ by terms which only 
become zero in the N 00 limit. 



VI. BOSON-FERMION SYMMETRY AND FOUR-POINT Om CORRELATORS IN QH SYSTEMS 

In this section, we will consider a planar system of electrons in the presence of a magnetic field B and impurities 
represented by a random potential V{x). The corresponding Hamiltonian operator is, 

H = Ho + V , Ho=^{p-eAf, (6.1) 

where A is the vector potential (wc will use the symmetric gauge). 

Here, we can also introduce the quadratic operator according to eq. (2.4). As discussed, because of Wick's 
theorem, the four-point correlators associated with Om will encode important information about the mixed retarded 
and advanced products appearing in the density of localized states in eq. (1.2). In a similar way to eq. (2.12), we can 
introduce a single superfield with components (/?(x), V'(x) whose action is given by, 

5m = ^y (fx{^Om^ + iiOm^), (6.2) 

thus leading to a properly normalized representation of the four-point correlators containing Gm Green's functions. 

For a very strong magnetic field, the transitions between Landau levels are suppressed. Then, if we are interested 
in the case where the first Landau level is filled, we can consider the projection, 



(p = e ^^'^^''''u{z) , ip = e ^*'^^"v{z), 



(6.3) 
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FIG. 1: In these two figures the bold hne represents the correlator (5.12), for A'^ = 10 and = 70 respectively. The dashed 
line corresponds to the 1 — (sin^ s)/s^ law in eq. (5.11). We have considered 7 = 1 in the two figures 



where z = x + iy and k = eB/U. 

The action Sm for the projected fields can be written as, 

a{z,z) , a{z,z) = e~^'^^''''{uu + vv). (6.4) 

Wc would like to underline that eq. (6.4) is local in the random potential (containing up to quadratic terms). Then, 
if the potentials are uncorrelated at different points, the corresponding average will lead to an effective action which 
is local in the field a{z, z), 

S'J^ = \j dzdzgia), (6.5) 

where g{a) depends on the type of disorder. 

Then, we will be able to follow refs. [14, 15], defining a superfield <i>(z, 9) = u{z) + k9v{z) to write, 

a{z,z) = rf^d^e-5'''(^2+««)$i, (6.6) 
(the normalization for the Grassmann measure is / dOdS e~^^ = I/tt) and construct a supersymmetric representation, 

S^J^ = ^§ / dzdzdedOe-'^'^'^^-^+'^'^hm , h{a) = £ ^g{P), (6.7) 

with the remarkable boson- fermion symmetry appearing in ref. [15], when computing the total density of states in 
QH systems. This is an interesting rcsiilt as this is the underlying symmetry behind the exact expression for the total 
density. 

As discussed in that reference, if the traditional approach to treat mixed retarded and advanced correlators is 
followed, the boson-fermion symmetry is lost because of the superfield doubling used to represent each prescription. 
Here, we have seen that this important symmetry is restored if the single superfield representation for mixed correlators 
is instead considered. 



'J 



dzdz 



V ■ 



heB 
Am 



E 



(r/2 + ivif 



VII. CONCLUSIONS AND PERSPECTIVES 



In this work we have introduced a new representation for mixed retarded and advanced four-point correlators. 
These quantities are the relevant ones when studying localization and transport properties in disordered systems. On 
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FIG. 2: In this set of figures the bold line represents the correlator (5.12), for N = 10, 30, 50 and 70. The dashed Une corresponds 

to the 1 — (sin'^ law in cq. (5.11), while the dotted line corresponds to the N DOS-DOS correlator obtained in terms of 

orthogonal polynomials[18]. We have considered 7 = 1 in the four figures. 



the other hand, these are precisely the kind of correlators which are more difficult to deal with. This comes about 
from the mixed prescription involved, which leads to a doubhng of the field copies or the superfield representation in 
the traditional approaches. 

Here, we have considered a single Green's function Gm for an operator Om, quadratic in the system's Hamiltonian, 

which already contains the abovcmcntioncd mixing. In this manner, wc have been able to encode the nontrivial 
information in four-point correlators essentially using a single prescription. In general, this Om correlator differs from 
the mixed four-point correlator or the density- density correlator by only retarded and only advanced terms, which 
can be computed by standard procedures. 

These ideas have been tested in the simpler context of Random Matrix theory. For this aim, we have obtained in 
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section V, in the N ^ oo limit, a direct relationship between the density- density correlator, 

{p{E,)p{E2))/{p{E,)p{E2)), 

and our four-point correlator for Om operators, whose GUE average has been computed in section III. Our closed 
expression (5.12) for the (finite A'') 0„ correlator clearly shows the onset of the well known asymptotic level spacing 
correlation. 

Then, we have seen that the single superfield method is an interesting alternative to be considered when studying 
disorder. In particular, in the context of quantum Hall systems, we have shown in section VI that the boson-fermion 
symmetry is restored for correlators encoding information about mixed retarded and advanced Green's functions, thus 
opening the possibility of improving our understanding of localization [19]. 
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Appendix 

Here we show how to compute traces of matrices which are functions of K = <f (^Tp — tp 'ip. For a function f{x), 
with a corresponding expansion f{x) = J^nf^x"', we have, trf{K) = J^nf^trK"'. On the other hand, a direct 
computation shows that. 



trK = a-(3 
trK'^ = + 2^71 - 0- 
trK^ = + 3{a + f3)n]I- P^, (7.1) 



where, 

a = ip-Tp , f3 = ip-tp , n-ip-ip. (7.2) 
It is easy to see the general form of the trace for a given power of K, 

trK"" = a" + n(a"-2 + a""^/? + a"-*/?^ + . . . + a^/S""^ + a(3"-^ + P"-'^)iJ4l - /?", 
or in an equivalent form, 

trK"=a"+n ^ A*/!-/?" , n ^ 0. (7.3) 

a — p 

Then, we arrive to a formula that gives the trace of a function of K just in terms of the f7(A'')-invariants (7.2): 

trf{K) = tr (^foI + Y.JnK^^ 

n — 1 



= Nfo + f{a) - f{l3) + lli^)-JM^,-p. (7.4) 

a — p 
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